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We present an improved protocol for entanglement purification of bipartite mixed states. 
The protocol requires two-way classical communication and hence implies an improved 
lower bound on the quantum capacity with two-way classical communication of the 
quantum depolarizing channel. 



1 Introduction 

Quantum information theory and quantum computation study the use of quantum physics in 
information processing and computationflj. Many important results such as quantum tele- 
portation, superdense coding, factoring and search algorithms make use of quantum entangle- 
ments as fundamental resources 

0131 HE]. Therefore, entanglement purification protocols, 
the procedures by which we extract pure-state entanglements from mixed states, merit our 
study. 

In this work, we follow the framework of [6l [7] and present a new purification protocol 
with improved yields 

1 . 1 Notation 

We denote von Neumann entropy by S(p) and Shannon entropy by H(po,pi, . . .). The fol- 
lowing notation is used for the four Bell states: 
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To facilitate our discussion later on, we also use two classical bits to label each of the Bell 
states: 

$+ = 00 
= 01 
= 10 

* _ = 11 (2) 

and we concentrate on the generalization of the Werner state [5]: 

Pf = F |*+) <$+| + 1^ ( |$-) ($-| + |*+) (*+| + |*-) (*" | \ (3) 

This work concerns entanglement purification protocols. At the beginning of these protocols, 
two persons, Alice and Bob, share a large number of quantum states pp, say Pp n , and they 
are allowed to communicate classically, apply unitary transformations and perform projective 
measurements. We place no restriction on the size of their ancilla systems so that we lost 
no generality in restricting their local operations to unitaries and projective measurements. 
In the end, the quantum states T shared by Alice and Bob are to be a close approximation 
of the maximally entangled states (|3> + ) ($+|)® m , or more precisely we require the fidelity 
between T and ($ + |)® m approaches one as n goes to infinity. We define the yield of 

such protocols to be m/n. 

1.2 Previous protocols 

We briefly review some previous protocols in the following subsections: 
1.2.1 Universal hashing 

Universal hashing was introduced in [5) and requires only one-way classical communication. 
The scenario is the same as what was described in section 11.11 The hashing method works 
by having Alice and Bob each perform some local unitary operations on the corresponding 
members of the shared bipartite quantum states. They then locally measure some of the 
pairs to gain information about the Bell states of the remaining unmeasured pairs. It was 
shown that each measurement can be made to reveal almost 1 bit of information about the 
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unmeasured pairs; therefore, by measuring nS(pp) pairs, Alice and Bob can figure out the 
identities of the remaining unmeasured pairs with high probability. Once the identities of the 
Bell states are known, Alice and Bob can convert them into the standard states $ + easily. 
This protocol distills a yield Dh = [n — nS(pF)]/n = 1 — S(pf)- 



1.2.2 The recurrence method 
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Fig. 1. The recurrence method. 

The recurrence method[6j [9j is illustrated in figHJ Alice and Bob put the quantum states 
pf!" into groups of two and apply XOR operations to the corresponding members of the 
quantum states p® 2 , one as the source and one as the target. They then take projective 
measurements on the target states along the z-axis, and compare their measurement results 
with the side classical communication channel. If they get identical results, the source pair 
"passed" ; otherwise the source pair "failed" . Alice and Bob then collect all the "passed" pairs, 
and apply a unilateral 7r rotation a x followed by a bilateral tt/2 rotation £?J? This process is 
iterated until it becomes more beneficial to pass on to the universal hashing. If we denote the 
quantum states by p = p 00 |$ + ) ($ + | + Poi \^ + ) (* + | + pw |$~) (*~| + Pn I* - ) (*~|, then 
this protocol has the following recurrence relation: 



POO — {POO + Pw)/Ppass', Poi 

Pw = 2poiPn/P: 



pass ■) 



p'n 



(POI +P 2 ll)/Ppas S ; 
ZpOOPw/Ppass] 



and 



Ppas 



Poo + Poi + Pw + Pii + 2pooPio + 2p iPn 



c As mentioned in [6], the application of a a x and B x rather than a twirl was proposed by C. Macchiavello. 
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1.2.3 The Maneva-Smolin method 

The Maneva-Smolin method [7] is illustrated in figlU Alice and Bob first choose a block size to 
and put the quantum states into groups of to. They then apply bipartite XOR gates between 
each of the first m — 1 pairs and the mth pairs. After that, they take measurements on 
these mth pairs along the z-axis, and compare their results with side classical communication 
channel. If they get identical results, they perform universal hashing on the corresponding 
to — 1 remaining pairs; if they get different results, they throw away all to pairs. The yield 
for this method is: 




Bob 

Fig. 2. The Maneva-Smolin method when m = 4. 

2 Entanglement purification protocol 

In section 12.11 we will present our new entanglement purification protocol and compare its 
yield with the yields obtained by the recurrence method [6] and the Maneva-Smolin method 
[7J. In section [2T21 we will give a closed-form expression for the yield of this new protocol. 

2.1 New protocol and improved yield 

Our protocol is illustrated in figEl Alice and Bob share the quantum states pf, n and put 
them into groups of four. They then apply the quantum circuit shown in fig(3] and take 
measurements on the third and fourth pairs along the x- and z-axis respectively. Using the 
side classical communication channel, they can compare their results with each other. If they 
get identical results on both measurements, they keep the first and second pairs and apply 
universal hashing[6 . If either of the two results disagrees, they throw away all four pairs. 

The four pairs can be described by an 8-bit binary string, and since these are mixed 
states they are in fact probability distribution over all 256(= 2 8 ) possible 8-bit binary strings. 
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Fig. 3. Alice and Bob put the quantum states p|?" into groups of four, and apply quantum circuits 
consisting only of XOR gates. They then take measurements on the third pair along the x-axis and 
the fourth pair along the z-axis. With side classical communication, they compare their results and 
if both results agree, they apply universal hashing on the first two pairs. If either measurement 
result disagrees, they throw away all four pairs. The yield of this protocol is plotted on fig. [4] 

The quantum circuit consists only of XOR gates and therefore maps the 8-bit binary strings, 
along with their underlying probability distribution, bijectively to themselves. Let us call 
these probability distributions P(aici2&i£»2Ci 02^1^2) and P 1 (aia2bib 2 cic 2 did 2 ) ■ 

Our quantum measurements on the third and fourth pairs are simply looking at the 5th 
bit (measurement on the third pair along x-axis) and the 8th bit (measurement on the fourth 
pair along z-axis), where a "0" means Alice and Bob getting identical results and a "1" means 
their getting opposite results. For example, if the 8-bit binary string is "0102616201 C2<iid2 = 
00100111", which corresponds to the quantum states $ + $ _ \f ,+ \f ,_ , then Alice and Bob will 
get identical results on the third pair but opposite results on the fourth. The "pass" prob- 
ability is Pp asa = J2 ai , a 2 , 61, b 2 c 2 , die{o, 1} -P'(ai a 20i62Oc 2 diO) and the post-measurement proba- 
bility distribution is Q(aia 2 &i&2) = 12 C2 , die {0.1} P' ( a i a 2bib 2 0c2di0) /p pass - The yield of this 
method [7] is: 



Ppass A _ H(Q(aia 2 bib2)) 
2 \ 2 

where H ((5(01026162)) is the Shannon entropy function. Fig. [4]compares the yield of our new 
method with the recurrence method and the Maneva-Smoline method. 




2.2 Closed-form expression 

The quantum circuit that Alice and Bob apply to the quantum states p^ 4 consists only of 
XOR gates and therefore maps the 8-bit binary strings bijectively to themselves. We call this 
bijection /: 
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Recurrence Maneva-Smolin New method 




Initial Fidelity 



Fig. 4. The dotted line is the yield for modified recurrence method [6]; the dash line is for the 
Maneva-Smolin method [7]. The yield of our new method is represented by the solid line, and 
there is an improvement over the previous methods when the initial fidelity is between 7.5 and 
8.45. 
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/:{0,1} 8 — > {0,1} 8 
(ai,a 2 ,bi,b 2 ,ci,C2,di,d 2 ) i — ► (oi © di, a 2 © c 2 , 61 © di, 62 © c 2 , 

ai © 61 © ci © di,c 2 ,di,a 2 © & 2 © c 2 © d 2 ) 



In table [U we list the quantum states that lead to identical measurement results for Alice and 
Bob and their probabilities in the ensemble p® 4 . Therefore, we can write down expressions 
for the terms p paS s and i?((5(aia 2 6i6 2 ) in equation^ as follows: 



Ppass = F 4 + 18F 2 G 2 + 24FG 3 + 21G 4 (5) 

„,_, , , /F 4 + 3G 4 \. /F 4 + 3G 4 \ /2F 2 G 2 + 2G 4 \ , (2F 2 G 2 + 2G i 
H{Q{a 1 a 2 b 1 b2)) = -[ J log 2 ( J - 91 1 log 2 f 

^ Ppass ' ^ Ppass ' ^ Ppass ' ^ Ppass 

„/4FG 3 \ MFG 3 \ 



where G = (1 -F)/3. 
3 Conclusions 

We presented a new protocol for entanglement purification assisted by two-way classical com- 
munication. It was shown in [6] that such a protocol corresponds to quantum capacity assisted 
by two-way classical communication of the quantum depolarizing channel, and hence we have 
a new lower bound for this capacity. In section 12.11 we applied the new protocol to the 
quantum states 



PF = F\$+) ($+| + ^— ^( |$-> ($-| - - hi' ; ;* 



however, our method also works for any Bell-diagonal quantum states 



P = Poo |$ + > ($ + | +P01 |* + ) (* + | + Pio I* - ) (*"| + Pll |*~) (*"| 



Equation ([5]) and ([6]) then become 
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Ppass = (Poo + POI + PlO + Pll) + 6 X 4 PooPOlPlOPll + 3 X ^ 2 PiP. 



ije{o,i} 2 



'Poo + Pm + Pio + Pn \ /Poo + Poi + Pio + Pii 



H(Q( ai a 2blb2 )) = ( Poo+Poi+Pio+Pn \ lQg2 / 



-6 x 



-3 x 



Ppa ss ' ^ Ppa s s 

4pooPoiPioPii \ / 4poopoipiopii \ 

Ppass ' ^ Ppass ' 

2 PooPoi + 2 P?oP?i \ , / 2pg pg 1 + 2p? pf j 



/ ^PooP5i +^PioPn \ log / 

^ T)nnss ' ^ 



Ppass ' ^ Ppas 



_ 3 x / gPiopjo + 2pjip|i \ (IploPlo + 2poiP?r 



Ppass ' * Ppass 

3 x / gPogPu^jPoiPio N bg /2pg P?i + 2p^pf 



Ppass ' ^ Ppass 

With these equations, we can combine the recurrence method and our new method: we start 
with the recurrence method and pass on to our new method rather than universal hashing. 
Indeed, there are improvements, but they occur over segments of narrow regions and the 
improvements are insignificant. Therefore we believe these improvements have only to do 
with the number of recurrence steps performed before passing on to universal hashing, and 
we will spare the readers with the details. 

An obvious direction is to look for new protocols with better yields. In particular, as 
was raised by E.N. Maneva and J. A. Smolin in [7], is there a way to iterate this protocol 
rather than passing on immediately to universal hashing? It is conceivable that, with such 
an adaption, we can have improvements over a wider range of initial fidelity. 

After the completion of this work, it came to our attention similar works have been carried 

out in nnnnj. 
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Table 1. The quantum states that lead to identical results for Alice and Bob. 
G = (1 - F)/3; $+ = 00; = 01; *~ = 10; = 11; 



P \a,\CL2b\b2C\C2d\d'2) 




/ (&l 6t2 b\ b2 ci C2 d\ d,2 ) 


tr c (i (/ \a1a2b1b2c1c2d1d2)) 


F 4 


00000000 


00000000 


0000 


G 4 


01010101 


00000100 


0000 


G 4 


10101010 


00000010 


0000 


G 4 


11111111 


00000110 


0000 


F'^G-* 


00010001 


00010000 


0001 


F 2 G 2 


01000100 


00010100 


0001 


G 4 


10111011 


00010010 


0001 


G 4 


11101110 


00010110 


0001 


F Z G 2 


00101000 


00100000 


0010 


G 4 


01111101 


00100100 


0010 


F 2 G 2 


10000010 


00100010 


0010 


G 4 


11010111 


00100110 


0010 


FG 3 


00111001 


00110000 


0011 


FG 3 


01101100 


00110100 


0011 


FG 3 


10010011 


00110010 


0011 


FG 3 


11000110 


00110110 


0011 


F 2 G 2 


00010100 


01000100 


0100 


F 2 G 2 


01000001 


01000000 


0100 


G 4 


10111110 


01000110 


0100 


G 4 


11101011 


01000010 


0100 




00000101 


01010100 


0101 


F 2 G 2 


01010000 


01010000 


0101 


G 4 


10101111 


01010110 


0101 


G 4 


11111010 


01010010 


0101 


F 2 G^ 


00111100 


01100100 


0110 


G 4 


01101001 


01100000 


0110 


G 4 


10010110 


01100110 


0110 


F 2 G 2 


11000011 


01100010 


0110 


FG 3 


00101101 


01110100 


0111 


FG 3 


01111000 


01110000 


0111 


FG 3 


10000111 


01110110 


0111 


FG 3 


11010010 


01110010 


0111 


F 2 G 2 


00100010 


10000010 


1000 


G 4 


01110111 


10000110 


1000 


F 2 G 2 


10001000 


10000000 


1000 


G 4 


11011101 


10000100 


1000 


F 2 G 2 


00110011 


10010010 


1001 


G 4 


01100110 


10010110 


1001 


G 4 


10011001 


10010000 


1001 


F 2 G 2 


11001100 


10010100 


1001 


F^G* 


00001010 


10100010 


1010 


G 4 


01011111 


10100110 


1010 


F 2 G 2 


10100000 


10100000 


1010 


G 4 


11110101 


10100100 


1010 


FG 3 


00011011 


10110010 


1011 


FG 3 


01001110 


10110110 


1011 


FG 3 


10110001 


10110000 


1011 


FG 3 


11100100 


10110100 


1011 


FG 3 


00110110 


11000110 


1100 


FG 3 


01100011 


11000010 


1100 


FG 3 


10011100 


11000100 


1100 


FG 3 


11001001 


11000000 


1100 


FG a 


00100111 


11010110 


1101 


FG 3 


01110010 


11010010 


1101 


FG 3 


10001101 


11010100 


1101 


FG 3 


11011000 


11010000 


1101 


FG 3 


00011110 


11100110 


1110 


FG 3 


01001011 


11100010 


1110 


FG 3 


10110100 


11100100 


1110 


FG 3 


11100001 


11100000 


1110 


F' Z G 2 


00001111 


11110110 


1111 


G 4 


01011010 


11110010 


1111 


G 4 


10100101 


11110100 


1111 


F 2 G 2 


11110000 


11110000 


1111 



